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20. 

I.  Ab-initio  quantum  chemical  investigations  of  the  mechanism  of  cationic 

polymerization,  (commenced  October  1,  1979) 

II.  Development  of  theoretical  procedures  for  prediction  of  crystal  structure: 

and  densities,  (commenced  April  1,  1980) 

I.  In  the  research  on  part  I  of  this  project  we  are  carrying  out  ab-initio 
quantum  chemical  calculations  and  generating  electrostatic  molecular  potential 
maps  for  cationic  polymerization:  on  energetic  monomers,  on  their  interaction 
with  catalysts  and  on  the  attack  of  the  monomer-catalyst  complexes  on  the 
substrates.  This  research  is  in  support  of  an  ongoing  overall  project  of  the 
Power  Program  Branch  of  ONR  designed  to  develop  energetic  binders  for  use 
with  energetic  propellants.  The  propensity  of  polymerization  of  more  conven¬ 
tional  cyclic  ethers  had  been  shown  experimentally  to  be  proportional  to  their 
basicity  and  ring  strain.  We  derived  a  theoretical  index  which  should  be 
proportional  to  the  basicity  of  the  monomers:  namely,  from  the  ab-initio 
quantum  chemical  calculations  on  the  monomers  generating  the  electrostatic 
molecular  potential  maps  in  space  around  the  monomer.  The  more  negative  the  map 
contours  the  farther  out  in  space,  the  more  basic  should  be  the  monomer. 

The  ab-initio  calculations  were  calculated  with  our  own  programs  which 
also  incorporate  several  desirable  options  for  calculations  on  large  molecules 
and  the  maps  were  generated  from  these  ab-initio  wave  functions. 

The  electrostatic  molecular  potential  contour  maps  generated  around  the 
exotically  substituted  cyclic  ether  monomers  proved  to  indicate  correctly 
their  different  propensities  for  polymerization  when  the  monomers  were  subse¬ 
quently  synthesized  and  polymerized  experimentally  by  the  other  researchers  on 
this  overall  ONR  project.  These  electrostatic  molecular  potential  contour 
maps  can  also  shed  insight  on  the  question  raised  as  to  good  copolymer 
candidates. 

We  have  also  carried  out  some  necessary  program  development  activities  in 
connection  with  this  research  in  improving  convergence  stability  of  our  current 
ab-initio  programs  since  this  is  sometimes  a  problem  for  nitro,  fluoro  and 
azido  substituents  and  implementing  new  approaches  to  forming  and  solving  SCF 
equations  for  reasonably  large  and  closely  related  molecules  or  molecular 
systems. 

II.  We  are  pursuing  several  approaches  to  the  prediction  of  crystal  structures 
and  densities.  The  first  is  to  calculate  (from  calculations  on  the  super¬ 
molecules)  the  SCF  intermolecular  interaction,  separate  it  into  its  various 
components  and  use  these  to  construct  the  crystal  structure. 

The  second  approach  is  to  carry  over  the  concept  recently  used  success¬ 
fully  in  calculations  for  biological  systems  where  the  same  ab-initio  calcul¬ 
ations  between  two  or  several  of  such  molecules  for  a  large  variety  of  inter¬ 
molecular  configurations  is  used  to  derive  atom  class-atom  class  pair-pair  and 
three-body  potentials  which  are  then  used  to  calculate  the  intermolecular 
structure  of  an  arrangement  of  molecules. 
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QUANTUM  CHEMICAL  INVESTIGATIONS  OF  THE  MECHANISM  OF 
CATIONIC  POLYMERIZATION 


Joyce  J.  Kaufman,  Principal  Investigator 
CONCISE  SUMMARY 


I.  Ab-initio  quantum  chemical  investigations  of  the  mechanism  of 
cationic  polymerization,  (commenced  October  1,  1979) 

II.  Development  of  theoretical  procedures  for  prediction  of  crystal 
structures  and  densities,  (commenced  April  1,  1980) 


I.  In  the  research  on  part  I  of  this  project  we  are  carrying  out  ab-initio 
quantum  chemical  calculations  and  generating  electrostatic  molecular 
potential  maps  for  cationic  polymerization:  on  energetic  monomers,  on 
their  interaction  with  catalysts  and  on  the  attack  of  the  monomer-catalyst 
complexes  on  the  substrates.  This  research  is  in  support  of  an  ongoing 
overall  project  of  the  Power  Program  Branch  of  ONR  designed  to  develop 
energetic  binders  for  use  with  energetic  propellants.  The  propensity  of 
polymerization  of  more  conventional  cyclic  ethers  had  been  shown  experi¬ 
mentally  to  be  proportional  to  their  basicity  and  ring  strain.  V/e 
derived  a  theoretical  index  which  should  be  proportional  to  the  basicity 
of  the  monomers:  namely,  from  the  ab-initio  quantum  chemical  calculations 
on  the  monomers  generating  the  electrostatic  molecular  potential  maps  in 
space  around  the  monomer.  The  more  negative  the  map  contours  the  farther 
out  in  space,  the  more  basic  should  be  the  monomer. 

The  ab-initio  calculations  were  calculated  with  our  own  programs 
which  also  incorporate  several  desirable  options  for  calculations  on 
large  molecules  and  the  maps  were  generated  from  these  ab-initio  wave 
functions. 

The  electrostatic  molecular  potential  contour  maps  generated  around 
the  exotically  substituted  cyclic  ether  monomers  proved  to  indicate  cor¬ 
rectly  their  different  propensities  for  polymerization  when  the  monomers 
were  subsequently  synthesized  and  polymerized  experimentally  by  the  other 
researchers  on  this  overall  ONR  project.  These  electrostatic  molecular 
potential  contour  maps  can  also  shed  insight  on  the  question  raised  as  to 
good  copolymer  candidates. 

We  have  also  carried  out  some  necessary  program  development  activities 
in  connection  with  this  research  in  improving  convergence  stability  of 
our  current  ab-initio  programs  since  this  is  sometimes  a  problem  for  nitro, 
fluoro  and  azido  substituents  and  implementing  new  approaches  to  forming 
and  solving  SCF  equations  for  reasonably  large  and  closely  related 
molecules  or  molecular  systems. 
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II.  We  are  pursuing  several  approaches  to  the  prediction  of  crystal 
structures  and  densities.  The  first  is  to  calculate  (from  calculations 
on  the  supermolecules)  the  SCF  intermolecular  interaction,  separate  it 
into  its  various  components  and  use  these  to  construct  the  crystal 
structure. 

The  second  approach  is  to  carry  over  the  concept  recently  used 
successfully  in  calculations  for  biological  systems  where  the  same 
ab-initio  calculations  between  two  or  several  of  such  molecules  for  a 
large  variety  of  intermolecular  configurations  is  used  to  derive  atom 
class-atom  class  pair-pair  and  three-body  potentials  which  are  then 
used  to  calculate  the  intermolecular  structure  of  an  arrangement  of 
molecules. 


» 
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QUANTUM  CHEMICAL  INVESTIGATIONS  OF  THE  MECHANISMS  OF 
CATIONIC  POLYMERIZATION 


Joyce  J.  Kaufman,  Principal  Investigator 


SUMMARY 


I.  Ab-initio  quantum  chemical  investigations  of  the  mechanism  of 
cationic  polymerization,  (commenced  October  1,  1979). 

II.  Development  of  theoretical  procedures  for  prediction  of  crystal 
structures  and  densities,  (commenced  April  1,  1980). 


I.  In  the  research  on  Part  I  of  this  project  we  are  carrying  out  ab-initio 
quantum  chemical  calculations  and  generating  electrostatic  molecular 
potential  maps  for  cationic  polymerization:  on  energetic  monomers,  on 
their  interaction  with  catalysts  and  on  the  attack  of  the  monomer-catalyst 
complexes  on  the  substrates.  This  research  is  in  support  of  an  ongoing 
overall  project  of  the  Power  Program  Branch  of  ONR  designed  to  develop 
energetic  binders  for  use  with  energetic  propellants.  A  binder  which 
itself  was  energetic  would  contribute  significantly  to  the  propellant 
effecti vi ty/weiqht  ratio.  The  compounds  of  choice  were  cyclic  ethers 
(preferably  3-  and  4-  membered  rings)  substituted  with  exotic  energetic 
substituents  such  as  fluoro-,  nitro-,  azido-,  azi domethyl- ,  etc.,  either 
singly  or  multiply  substituted.  The  primary  initial  Question  asked  was 
could  one  predict  the  relative  propensities  of  cationic  polymerization 
for  these  exotic  monomers  --  even  before  the  monomers  were  synthesized. 

The  only  previous  old  experimental  Japanese  paper  reported  that  for  cyclic 
ethers  substituted  with  more  conventional  substituents  (alkvl-,  alkoxy-, 
aryl-)  the  rate  of  polymerization  was  proportional  to  the  basicities 
of  the  monomers  and  the  ring  strain.  We  derived  a  theoretical  index  which 
should  be  proportional  to  the  basicity  of  the  monomers:  namely,  from  the 
ab-initio  guantum  chemical  calculations  on  the  monomers  generating  the 
electrostatic  molecular  potential  maps  in  space  around  the  monomer.  The 
more  negative  the  map  contours  the  farther  out  in  space,  the  more  basic 
should  be  the  monomer. 

The  ab-initio  calculations  were  calculated  with  our  own  programs 
which  also  incorporate  several  desirable  options  for  calculations  on 
large  molecules:  ab-initio  effective  core  model  potentials  (MODPOT)  which 
enable  calculations  of  the  valence  electrons  only  explicitly,  yet 
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accurately;  an  integral  prescreening  evaluation  (which  we  named  VRDDO  -- 
variable  retention  of  diatomic  differential  overlap)  especially  effective 
for  spatially  extended  systems; and  a  MERGE  technique  which  allows  efficient 
reuse  of  integrals  from  common  skeletal  fragments.  These  electrostatic 
molecular  potential  contour  maps  were  generated  around  the  exotical ly 
substituted  monomers  and  some  of  the  more  conventionally  substituted  mon¬ 
omers  which  had  earlier  been  studied  experimentally. 

The  electrostatic  molecular  potential  contour  maps  generated  around 
the  exotical ly  substituted  cyclic  ether  monomers  proved  to  indicate  correctly 
their  different  propensities  for  polymerization  when  the  monomers  were 
subsequently  synthesized  and  polymerized  experimentally  by  the  other 
researchers  on  this  overall  ONR  project.  These  electrostatic  molecular 
potential  contour  maps  can  also  shed  insight  on  the  question  raised  as  to 
good  copolymer  candidates.  Since  the  rates  of  polymerization  are  governed 
by  basicity  which  in  turn  are  predictable  from  these  maps,  the  more 
closely  these  maps  match  for  two  different  monomers,  the  more  well  suited 
these  two  monomers  are  for  copolymer  partners. 

We  have  also  carried  out  some  necessary  program  development  activities 
in  connection  with  this  research.  These  have  encompassed  two  areas: 
improving  the  convergence  stability  of  our  current  ab-initio  MODPOT/VRDDO/MERGE 
program  systems,  since  convergence  stability  can  sometimes  be  a  problem 
particularly  for  non-equilibrium  orientations  of  molecules  containing 
fluoro,  nitro,  or  azido  groups,  and  implementing  new  approaches  to  forming 
and  solving  SCF  equations  for  reasonably  large  and  closely  related  molecules 
or  molecular  systems. 


See  next  page  for  summary  of  Part  II. 
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II.  The  arrangement  of  molecules  in  crystals  is  important  in  understanding 
intermolecular  interactions  which  occur  in  crystals  during  detonation 
phenomena.  Crystals  densities  are  one  of  the  important  factors  influencing 
detonation  pressure  and  detonation  velocity.  Crystal  structures  and 
crystal  densities  can  be  calculated  by  ab-initio  MODPOT/VRDDQ/MERGE  SCF 
supermolecule  calculations,  supplemented  by  an  estimate  of  the  dispersion 
energies  (which  are  known  to  be  much  smaller  than  the  other  terms  for 
nitro  compounds).  However,  it  would  be  most  desirable  to  be  able  to 
calculate  crystal  structures  and  crystals  densities  by  a  intermolecular 
perturbation  technique. 

We  are  pursuing  several  approaches.  The  first  approach  is  to  calculate 
the  SCF  intermolecular  interaction  between  two  or  several  such  molecules, 
then  to  separate  this  interaction  into  its  electrostatic,  polarization, 
exchange  and  charge-transfer  components  and  then  to  use  these  components 
(with  an  estimate  of  the  dispersion  energy)  to  construct  the  crystal 
structure.  We  have  done  the  analysis  for  this  approach  and  are  implementing 
it. 


The  second  approach  which  we  have  just  initiated  is  to  carry  over 
the  concept  recently  used  successfully  in  calculations  for  biological 
systems  where  the  same  ab-initio  calculations  between  two  or  several  of 
such  molecules  for  a  large  variety  of  intermolecular  configurations  is 
used  to  derive  atom  class-atom  class  pair-pair  and  three-body  potentials 
wnich  are  then  used  to  calculate  the  intermolecular  structure  of  an 
arrangement  of  molecules. 


I 

t 

i 
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I.  Ab-Initio  Quantum  Calculations  and  Electrostatic  Molecular  Potential 
Maps  For  the  Mechanism  of  Cationic  Polymerization 

By  making  energetic  binders  for  propellants,  the  efficiency/weight 
ratio  of  propellants  could  be  increased  significantly. 

ONR  Power  Program  Branch  is  sponsoring  an  overall  effort  in  this  area 
which  combines  other  groups  in  experimental  synthesis  and  polymerization 
with  our  theoretical  group.  The  objective  of  the  QNR  project  is  to  prepare 
energetic  monomers  of  cyclic  ethers  (of  ring  sizes  three-to-five-membered) 
substituted  with  energetic  groups  (such  as  nitro-,  fluoro-,  aziao-, 
azidomethyl-)  both  single  substituents  and  mixed  substituents.  We  were 
asked  by  ONR  what  insight  quantum  chemical  calculations  could  shed  on 

1)  the  question  of  the  propensity  of  the  cyclic  ethers  with 
energetic  substituents  to  undergo  cationic  polymerization  prior 
even  to  synthesis  of  the  energetic  monomers  themselves; 

2)  what  happens  as  a  catalyst  (a  Lewis  acid)  approaches  the  cyclic 
ether  monomer  and  what  might  be  preferable  catalysts; 

3)  what  happens  as  a  complex  of  cyclic  ether  +  catalyst  approaches 
another  molecule  of  cyclic  ether  monomer; 

4)  what  would  be  good  partners  for  co-polymerization  in  the  event 
the  polymer  from  a  single  monomer  was  too  refractory  --  too 
brittle. 

It  had  been  observed  by  Japanese  investigators  (Aoki,  et  al.),  who 
earlier  studied  experimentally  three-to-six  membered  cyclic  ethers  with 
more  conventional  substituents  (alkyl,  alkoxy,  aryl)  that  the  rate  of 
polymerization  of  cyclic  ethers  was  proportional  to 

a)  the  base  stre  gth  of  the  cyclic  ether 

b)  the  ring  strain  of  the  cyclic  ether. 

Thus,  to  determine  propensity  for  polymerization,  we  had  to  come  up 
with  theoretical  quantum  chemical  indices  that  would  be  predictive  of 
basicity  and  of  ring  strain.  Our  previous  ab-initio  quantum  chemical 
calculations  on  morphine  (an  N-methyl  narcotic  agonist)  and  on  nalorphine 
(an  N-allyl  narcotic  agonist-antagonist)  had  shown  that  charge  alone  on  an 
atom  does  not  seem  to  be  sufficiently  sensitive  to  distinguish  different 
degrees  of  basicity  especially  for  closely  related  molecules.  However, 
the  use  of  electrostatic  molecular  potential  contour  maps  did  hold  promise 
of  being  a  sensitive  criterion  for  basicity.  Comparisons  of  negative 
regions  of  the  maps  around  the  0  atom  in  the  cyclic  ether  ring  (both 
spatial  extent  and  magnitude)  should  be  indicative  of  the  relative  basicities 
of  these  exotic  monomers.  For  the  cyclic  ether  monomers,  we  reasoned 
that  the  larger  in  space  the  negative  contours  of  the  map  around  the  oxygen 
atom  in  three  dimensions,  the  more  basic  the  monomer  would  be  and  the 
more  prone  to  cationic  polymerization.  For  each  monomer  we  are  carrying 
out  the  ab-initio  calculations  to  optimize  the  geometry  of  the  monomer. 

As  the  geometry  is  optimized  we  then  generate  electrostatic  molecular 
potential  contour  maps  in  three  dimensions  around  the  monomer  from  the 
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calculated  electronic  wave  functions  of  the  monomer.  These  maps  indicate 
vividly  the  basicity  of  the  oxygen,  from  which  direction  a  cation  will 
prefer  to  attack  the  monomer  and  are  especially  valuable  in  comparing 
three  dimensional  relative  basicities  of  various  monomers.  These  maps 
have  been  generated  around  a  number  of  the  monomers  and  around  some 
monomer-catalyst  complexes. 

We  are  carrying  out  ab-initio  MODPOT/VRDDO  and  ab-initio  MODPOT/VRDDO/MERGE 
calculations  for  3-,  4-,  5-  and  6-member  cyclic  ethers  with  various 
substituents.  The  3-  and  4-member  rings  are  the  ones  of  choice  for  the 
experimental  program.  In  particular,  the  four  member  rings  3,3-fluoro- 
nitrooxetane  and  3,3-dinitrooxetane  were  of  the  first  interest  experimentally. 

We  carried  out  ab-initio  MODPOT/VRDDO/MERGE  calculations  on  oxetane  and 
on  the  3,3-fluoronitrooxetane  and  3,3-dinitrooxetane  and  generated  electro¬ 
static  molecular  potential  maps  around  the  molecules  in  three  directions 
i)  parallel  to  the  mean  molecular  plane  (xz);  ii)  parallel  to  the  (xy) 
plane  and  iii)  parallel  to  the  (yz)  plane.  The  results  for 

(a)  Oxetane  (Figures  la  and  lb) 

(b)  3,3-Fluoroni trooxetane  (Figures  2 a  and  2b) 

(c)  3, 3-Dini trooxetane  (Figures  3a  and  3b) 

indicate  quite  clearly  that  oxetane  is  the  most  basic,  3,3-fluoronitrooxetano 
is  next  and  3,3-dinitrooxetane  is  the  least  basic.  The  ring  strains 
are  not  expected  to  differ  nearly  as  drastically  in  these  molecules  as  do 
ttie  basicities.  The  maps  suggested  that  the  3,3-dinitrooxetane  should 
be  so  weakly  basic  that  it  would  have  trouble  undergoing  cationic  polymer¬ 
ization.  Very  recent  experimental  data  which  indicates  it  appears  to 
be  very  difficult  to  polymerize  dinitrooxetane  validates  our  hypothesis. 

Thus,  it  now  seems  possible  that  by  carrying  out  the  ab-initio  quantum 
chemical  calculations  on  these  cyclic  ether  monomers  with  combinations 
of  exotic  energetic  substituents  even  prior  to  the  synthesis  of  the  monomers 
themselves,  we  will  be  able  to  get  considerable  insight  into  the  propensity 
of  these  monomers  to  undergo  cationic  polymerization.  We  are  currently 
developing  theoretical  procedures  for  calculating  relative  ring  strains 
in  these  cyclic  ethers.  Thus,  in  response  to  question  1),  it  would  be 
expected  that  oxetane  would  polymerize  more  easily  than  fluoroni trooxetane 
and  that  more  easily  than  dinitrooxetane. 

These  electrostatic  molecular  potential  contour  maps  can  also  be 
used  to  help  answer  question  4)  as  to  good  candidate  partners  for  copolymer¬ 
ization.  We  hypothesized  that  since  the  rate  of  polymerization  is  very 
dependent  on  the  basicity  of  the  cyclic  ethers,  for  good  copolymerization 
the  copolymer  molecules  should  have  reasonably  similar  basicities  (as  well 
as  ring  strains),  because  if  one  compound  was  much  more  basic  than  the 
other,  it  would  sel f-polymerize  so  fast  that  it  would  not  have  much 
tendency  to  copolymerize.  3,3-bis(azidomethyl Joxetane  forms  a  polymer 
which  is  too  refractory  (too  brittle).  The  question  was  asked  as  to 
possible  good  copolymer  candidates  for  3,3-bis(azidomethyl )oxetane. 

The  maps  showed  that  3,3-bis (azidomethyl joxetane  was  less  basic  than 
oxetane,  more  basic  than  3,3-fluoronitrooxetane  and  apparently  fairly 
similar  in  basicity  to  3-azidooxetane.  Thus,  3-azi dooxetane  might  be  a 
good  choice  for  a  copolymer  candidate. 
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To  shed  insight  on  questions  2)  and  3)  we  are  also  carrying  out 
ab-initio  M0DP0T/VRD00  and  ab-initio  MODPOT/VRDDO/MERGE  calculations  on 
BF3  approaching  various  series  of  substituted  cyclic  ethers  and  generating 
the  electrostatic  molecular  potential  contour  maps.  These  results  indicate 
the  charge  redistribution  and  changes  in  the  maps  as+the  reactants  approach. 
Similar  calculations  are  carried  out  for  attack  of  H  on  the  substituted 
cyclic  ethers  since  there  has  been  speculation  by  experimental  polymer 
chemists  that  perhaps  it  is  the  attack  of  H+  (formed  from  a  Lewis  acid-H  0 
complex)  that  triggers  the  polymerization.  However,  experimentally  it  has 
been  observed  that  different  Lewis  acid  catalysts  have  different  effects 
on  polymerization  of  the  same  monomer.  For  example,  it  is  possible  to 
polymerize  certain  of  these  energetic  monomers  with  PF5  but  not  with  BF3. 

This  would  seem  to  indicate  that  a  catalyst-monomer  complex  might  still 
be  important.  We  suggested  several  experiments  that  could  help  clarify 
the  picture.  First,  to  measure  the  differing  amounts  of  H+  formed  by  various 
Lewis-acid  catalysts  in  the  polymerization  solvent  mixture  (without  any 
monomer  as  well  as  with  monomers)  to  see  if  the  amount  of  lr  formed  could 
explain  the  difference  in  propensity  of  a  monomer  to  polymerize  with+ 
different  Lewis-acid  catalysts.  Second,  to  add  various  amounts  of  H  to 
solutions  of  monomer  plus  ineffective  Lewis-acid  catalyst  to  see  if  the  + 
extra  H+  could  trigger  the  reaction.  Third,  to  add  various  amounts  of  H 
to  solutions  of  monomer  plus  effective  Lewis-acid  catalysts  to  see  if  the 
rate  or  character  of  the  polymerizations  are  affected. 


,*t& 1 


UiW-'.'-v***  <e«*  f  -  a* 


Molecular  Calculation  Summary 


The  parent  molecules  studied  so  far  have  come  from  three  sources: 

•  ONR  molecules  of  interest 

•  Molecules  from  the  Aoki,  Harita  et  al.  reaction  rate  study 
[J.  Polymer  Sci.  Al  6,  2585  (1968)] 

•  A  systematic  study  of  single  alpha  substituents  to  round 
out  and  extend  the  Aoki  et  al .  results. 

These  parent  molecules  are  also  being  investigated  in  attack  by  a 
proton  and  by  boron  trifluoride. 

The  following  abbreviations  are  used  in  the  table: 

for  origin:  N  =  ONR  molecule  of  interest. 

J  =  Aoki  et  al.  study. 

S  =  Systematic  single  a  substituent  study. 


for  status:  0  =  Geometry  optimization  completed. 

A  =  Geometry  optimization  approaching  completion. 

G  =  Geometry  optimization  being  done. 

C  =  Convergence  difficulties  encountered,  but 
stability  problems  are  in  the  process  of 
being  overcome. 

I(x)  =  In  plane  proton  attack  (where  the  reacting 

species  approaches  in  the  plane  of  the  parent 
that  contains  the  hetero-oxygen  and  its  ring 
neighbors) . 

Geometry  optimization  completed,  x  is  the 
equilibrium  H+  -  0  distance. 

IA(x)  =  Similar  to  I (x)  but  geometry  optimization 
approaching  completion. 

PA(x)  =  Out  of  plane  proton  attack  geometry  optimization 
x  is  the  angle  of  optimal  attack. 
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Calculational  Status 

El ectrostatic  Molecular  Potential  Contour  Maps  Calculated 

Oxirane  with  BF3  distances  of  1.0,  1.5,  2.0,  2.5,  5.0  A 

2- Hydroxyoxirane 
Oxetane 

3. 3- Di ni trooxetane 

3-  FI  uoro-3-Ni  trooxetane 
3-Azidooxetane 

3-Azi dooxetane  with  BF3  B-0  distances  of  1.0,  2.0,  2.5,  5.0  A 

3.3- Bis(azidomethyl Joxetane 
Oxolane  (THF) 

Oxolane  with  BF3  B-0  distance  of  1.0  A 
Oxane 

Oxane  with  BF3  B-0  distances  of  1.5,  2.0,  5.0  A 
4,7-Bis(fluorodinitromethyl )dioxepane 
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We  have  also  carried  out  some  necessary  program  development  activities 
in  connection  with  this  research.  These  have  encompassed  two  areas. 

We  have  improved  the  convergence  stability  of  our  ab-initio  MODPOT/ 
VRDDO/MERGE  program  systems  since  SCF  convergence  stability  can  sometimes 
be  a  problem,  particularly  for  non-equilibrium  orientations  of  molecules 
containing  fluoro,  azido  or  nitro  substituents.  A  simple,  but  elegant, 
attempt  at  stability  control  is  to  monitor  the  quality  of  the  density 
matrix.  Normally  the  reduced  density  matrix  is  computed  as  a  sum  over  all 
occupied  orbitals.  In  a  destabilized  iteration,  a  normally  occupied  orbital 
can  become  unoccupied.  In  our  new  approach,  this  orbital  is  not  included 
in  the  density  sum  and  the  density  is  weighted  to  correct  for  its  absence. 
This  prevents  the  contamination  of  unoccupied  orbitals  and  contributes  to 
iterative  stability.  In  essence,  this  is  a  ersatz  projection  technique 
and  works  well  for  these  systems. 

The  improvement  to  traditional  SCF  methodology  is  based  on  the 
chemical  knowledge  that  if  one  part  of  the  environment  of  a  reasonably 
large  molecule  changes,  then  radical  restructuring  of  the  electronic 
structure  will  most  likely  be  confined  to  the  locality  of  the  change.  This 
is  reflected  in  changes  to  the  coefficients  of  atomic  basis  functions  in 
the  vicinity  of  the  site.  We  have  formulated  an  exact,  variational  pro¬ 
cess  for  freezing  in  the  molecular  wave  function  the  coefficients  of  some 
of  the  atomic  basis  functions  (derived  as  a  result  of  a  previous  molecular 
calculation)  and  iteratively  computing  the  unfrozen  "reactions"  functions 
coefficients.  In  practice,  however,  one  needs  the  entire  set  of  coefficients 
unfrozen,  so  we  have  devised  and  implemented  a  revised  SCF  scheme  whereby 
some  coefficients  are  frozen,  SCF  convergence  achieved  for  the  unfrozen  set, 
and  then  converging  the  entire  unfrozen  set.  The  potentialities  for  com¬ 
putational  enhancement  are  great  since  the  effective  basis  set  size  is 
reduced  substantially  during  freezing.  Consider  a  calculation  of  BF3 
attacking  THF  using  a  small  basis  of  4  functions  per  atom,  except  for  hydro¬ 
gen,  which  has  one.  Freezing  all  functions  except  those  on  oxygen,  boron 
or  fluorine  leaves  20  unfrozen.  The  work  to  construct  the  Fock  matrix  to 
perform  a  diagonal ization  is  proportional  to  the  cube  of  the  size,  so  by 
freezing  we  expect  to  perform  only  9.4%  of  the  work.  If  we  additionally 
freeze  the  F  basis  functions,  the  work  is  now  only  0.6%  of  the  original! 

And  we  can  unfreeze  in  stages.  First  keep  unfrozen  B  and  0.  Then  unfreeze 
F,  and  finally  the  whole  system.  This  leads  not  only  to  speed  enhancement 
but  also  to  convergence  stability.  This  technique  will  be  particularly 
valuable  as  we  sequentially  bring  catalyst  up  to  monomer,  then  allow  the 
monomer-catalyst  to  attach  the  substrate  building  up  polymer  chains. 
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II.  Prediction  of  Crystal  Structure 

We  are  using  two  different  approaches  to  tackle  the  problem  of 
prediction  of  crystal  structure. 

In  the  first  approach,  while  crystal  structure  predictions  could 
be  done  by  carrying  out  ab-initio  SCF  calculations  of  the  entire  super¬ 
molecule,  it  is  more  desirable  to  try  to  carry  out  such  calculations  in 
terms  of  the  interactions  between  molecules  and  perturbation  theory. 

For  this  purpose  we  need  to  be  able  to  calculate  separately  the  electro¬ 
static,  exchange  (repulsion),  polarization  and  charge  transfer  energy 
contributions  to  intermolecular  interactions  (these  come  from  proper 
partitioning  of  the  SCF  solution)  and  to  estimate  accurately  the  dispersion 
energies,  and  we  are  working  on  these  aspects.  In  particular  the  polar¬ 
ization  term  must  be  modified  to  take  neighbor  effects  into  account. 

We  are  carrying  out  the  preliminary  analysis  and  are  coding  and  putting 
together  various  program  pieces  necessary  for  this  approach  for  ab-initio 
calculations  for  prediction  of  crystal  structures. 

Concomitantly,  through  the  graciousness  of  Dr.  P.  Claverie  (who 
spent  time  with  us  at  Johns  Hopkins  twice  last  year),  we  are  putting 
up  on  the  computers  we  use  a  less  rigorous  program  for  analysis  of 
crystal  geometries  and  for  intermolecular  interaction  energy.  This  was 
the  method  that  had  success  in  his  previous  analysis  of  the  crystal 
structure  of  nitrobenzene.  We  shall  explore  the  features  of  this 
method  while  we  are  writing  our  ab-initio  method. 

The  second  approach  which  we  have  just  initiated  is  to  carry  over 
the  concept  used  successfully  recently  in  biological  systems  of  carrying 
out  ab-initio  SCF  calculations  between  several  such  molecules  (in  a  variety 
of  orientations),  from  these  deriving  ab-initio  atom  class-atom  class 
pair-pair  and  three  body  potentials  which  can  then  be  used  to  calculate 
intermolecular  molecular  structure  accurately.  We  have  examined  these 
available  biological  ab-initio  potentials  and  they  lack  the  potentials 
for  nitro,  fluoro,  azido  and  for  aromatics  and  for  aliphatic  heterocyclics 
substituted  with  these.  Thus  we  shall  carry  out  the  ab-initio  inter¬ 
molecular  calculations  on  such  exotic  species  and  derive  the  necessary 
potentials.  A  large  number  of  such  ab-initio  calculations  are  necessary 
but  we  have  the  advantage  of  being  able  to  use  our  ab-initio  MODPOT/VRDDO/ 
MERGE  technique,  which  saves  considerable  time.  For  example,  for  nitro¬ 
benzene  which  was  our  original  test  molecule  the  ab-initio  calculation 
with  the  entire  6-3  G  basis  set  took  (on  a  CDC  6600  NOS/BE  1.3)  2691  seconds, 
the  VRD00  calculation  took  1000  secs.,  the  M0DP0T  calculation  took  717  secs., 
and  the  MODPOT/VRDDO  calculation  took  only  321  secs.  The  savings  in  time 
with  our  method  go  up  as  nx  as  the  molecules  get  larger  or  as  calculations  are 
carried  out  on  intermolecular  systems  (where  l<x<4  depending  on  the  size 
and  shape  of  the  molecule,  since  n4  is  the  ultimate  increase  in  the 
number  of  integrals  with  the  number  of  atomic  basis  functions  provided 
no  integrals  are  small).  Our  previous  extensive  testing  on  a  number  of 
molecules  had  shown  that  the  orbital  energies  and  gross  atomic  populations 
for  all  four  calculations  agreed  to  almost  0.001-0.002  a.u.  or  e,  and 
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potential  energy  curves  agreed  even  better,  to  0.001  -  0.0001  a.u.  As  an 
extra  bonus,  we  showed  last  year  that  the  M0DP0T  basis  sets  we  use 
are  so  well  balanced  that  there  is  very  little  superposition  error 
when  these  are  used  in  intermolecular  calculations. 

The  derived  atom  class-atom  pair  potentials  are  then  general  for 
a  particular  substituent  or  skeleton.  Thus,  they  can  then  be  used 
to  predict  crystal  structures  of  any  energetic  compound  containing 
those  structural  subunits. 
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III.  Joyce  J.  Kaufman  -  Invited  Lectures  Presented  on  this  ONR  Research 

A.  Invited  Lectures  Presented  at  the  Following  National  and 

International  Meetings,  Including  Results  on  the  Energetic  Monomers. 

March  1980.  International  Sanibel  Symposium  on  Atomic,  Molecular  and 
Solid  State  Theory,  Palm  Coast,  Florida. 

March  1980.  Symposium  on  Electrostatic  Potentials  in  Chemistry,  National 
American  Chemical  Society  Meeting,  Houston,  Texas.  In  press  in  the 
Proceedings  of  the  Meeting. 

June  1,  1980.  7th  Canadian  Theoretical  Chemistry  Conference,  Banff, 

A1 berta. 

July  1980.  NATO  Advanced  Study  Institute  on  Fast  Reactions  in  Energetic 
Systems,  Preveza  Beach,  Greece.  Also  presented  tutorial  lectures. 

In  press  in  the  Proceedings  of  the  Meeting. 


B.  Invited  Lectures  Presented  at  ONR  or  DOD  Meetings. 

August  1979.  Quantum  Chemical  Calculations  for  Energetic  Compounds 

(Explosives,  Propellants,  Etc.),  Conference  on  Thermal  Decomposition 
of  Propellants  and  Explosives,  Air  Force  Academy,  Boulder,  Colorado. 

September  1979.  Quantum  Chemical  Investigations  on  the  Mechanisms  of 

Cationic  Polymerization  and  Cationic  Graft  Polymerization;  Workshop 
on  Energetic  Binders,  ONR,  Washington,  D.C. 

April  1980.  Ab-Initio  Quantum  Chemical  Calculations  and  Electrostatic 
Molecular  Potential  Contour  Maps  Around  Energetic  Cyclic  Ethers. 

ONR  Contractors'  Meeting,  SRI,  Palo  Alto,  California. 

April  1980.  Ab-Initio  Quantum  Chemical  Techniques  for  Large  Molecules 
and  Application  to  Cationic  Polymerization.  Discussion  of  Navy 
research  interests.  Naval  Research  Laboratories. 


C.  In  addition,  at  the  request  of  Dr.  Richard  Miller,  our  ONR 
contract  monitor,  we  have  made  a  number  of  trips  to  ONR  and 
NRL  to  discuss  various  aspects  of  research. 
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IV.  Project  Personnel 

Joyce  J.  Kaufman,  Ph.D. ,  Principal  Investigator 
Frank  C.  Tobin,  Ph.D. 

The  majority  of  the  quantum  chemical  calculations  and  the  new 
SCF  program  development  reported  here  are  due  to  him. 

P.  C.  Hariharan,  Ph.D. 

Working  part  time  with  this  project  initially  since  he  was  already 
familiar  with  our  various  computer  programs  necessary  for  the  ab-initio 
quantum  chemical  calculations  and  electrostatic  molecular  potential 
contour  maps  and  he  has  also  made  various  addition  to  these  programs 
as  needed  for  the  present  research.  Now  participating  in  the  development 
of  the  program  for  predicting  crystal  structure. 
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University  of  California 
Livermore,  CA  94550 
Attn:  Dr.  R.  McGuire 

Lockheed  Missiles  and  Space  Co. 
P.0.  Box  504 
Sunnyvale,  CA  94088 
Attn:  Dr.  Jack  Linsk 
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Lockheed  Missile  &  Space  Co.  1 

3251  Hanover  Street 
Palo  Alto,  CA  94304 
Attn:  Dr.  H.  P.  Marshall 
Dept.  52-35 

Los  Alamos  Scientific  Lab  1 

P.0.  Box  1663 

Los  Alamos,  NM  87545 

Attn:  Dr.  R.  Rogers,  WX-2 

Los  Alamos  Scientific  Lab  1 

P.0.  Box  1663 

Los  Alamos,  NM  87545 

Attn:  Dr.  B.  Craig,  M  Division 

Naval  Air  Systems  Command  1 

Code  330 

Washington,  DC  20360 
Attn:  Mr.  R.  Heitkotter 
Mr.  R.  Brown 

Naval  Air  Systems  Command  1 

Code  310 

Washington,  DC  20360 
Attn:  Dr.  H.  Mueller 

Dr.  H.  Rosenwasser 

Naval  Explosive  Ordnance  1 

Disposal  Facility 
Indian  Head,  MD  20640 
Attn:  Lionel  Dickinson 
Code  D 

Naval  Ordnance  Station  1 

Code  5034 

Indian  Head,  MD  20640 
Attn:  Mr.  S.  Mitchell 

Naval  Ordnance  Station  1 

Code  PM4 

Indian  Head,  MD  20640 
Attn:  Mr.  C.  L.  Adams 

Dean  of  Research  1 

Naval  Postgraduate  School 
Monterey,  CA  93940 
Attn:  Dr.  William  Tolies 

Naval  Research  Lab  1 

Code  6510 

Washington,  DC  20375 
Attn:  Dr.  J.  Schnur 
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Naval  Research  Lab  1 

Code  6100 

Washington,  DC  20375 

Naval  Sea  Systems  Command  1 

Washington,  DC  20362 
Attn:'  Mr.  G.  Edwards,  Code  62R3 
Mr.  0.  Murrin,  Code  62R2 
Mr.  W.  Blaine,  Code  62R 

Naval  Sea  Systems  Command  1 

Washington,  DC  20362 
Attn:  Mr.  R.  Beauregard 
SEA  64E 


Naval  Surface  Weapons  Center  1 

Code  Rll 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  H.  G.  Adolph 

Naval  Surface  Weapons  Center  1 

Code  R13 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  R.  Bernecker 

Naval  Surface  Weapons  Center  1 

Code  R10 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  S.  J.  Jacobs 

Naval  Surface  Weapons  Center  1 

Code  Rll 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  M.  J.  Kamlet 

Naval  Surface  Weapons  Center  1 

Code  R04 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  D.  J.  Pastine 

Naval  Surface  Weapons  Center  1 

Code  R13 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Dr.  E.  Zimet 
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Naval  Surface  Weapons  Center  1 

Code  R101 

Indian  Head,  MD  20640 
Attn:  Mr.  6.  L.  Mackenzie 

Naval  Surface  Weapons  Center  1 

Code  R17 

Indian  Head,  MD  20640 
Attn:  Dr.  H.  Haiss 

Naval  Surface  Weapons  Center  1 

Code  Rll 

White  Oak,  Silver  Spring,  M0 
20910 

Attn:  Dr.  K.  F.  Mueller 

Naval  Surface  Weapons  Center  1 

Code  R16 

Indian  Head,  MD  20640 
Attn:  Dr.  T.  D.  Austin 

Naval  Surface  Weapons  Center  1 

Code  R122 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Mr.  L.  Roslund 

Naval  Surface  Weapons  Center  1 

Code  R121 

White  Oak,  Silver  Spring,  MD 
20910 

Attn:  Mr.  M.  Stosz 

Naval  Weapons  Center  1 

Code  3853 

China  Lake,  CA  93555 
Attn:  Dr.  R.  Atkins 

Naval  Weapons  Center  1 

Code  3205 

China  Lake,  CA  93555 
Attn:  Dr.  L.  Smith 

Naval  Weapons  Center  1 

Code  3205 

China  Lake,  CA  93555 
Attn:  Dr.  C.  Thelen 

Naval  Weapons  Center  1 

Code  385 

China  Lake,  CA  93555 
Attn:  Dr.  A.  Amster 


No.  Copies 

Naval  Weapons  Center  1 

Code  388 

China  Lake,  CA  93555 
Attn:  D.  R.  Derr 

Naval  Weapons  Center  1 

Code  388 

China  Lake,  CA  93555 
Attn:  Dr.  R.  Reed  Jr. 

Naval  Weapons  Center  1 

Code  385 

China  Lake,  CA  93555 
Attn:  Dr.  A.  Nielsen 

Naval  Weapons  Center  1 

Code  3858 

China  Lake,  CA  93555 
Attn:  Mr.  E.  Martin 

Naval  Weapons  Center  1 

China  Lake,  CA  93555 
Attn:  Mr.  R.  McCarten 

Naval  Weapons  Support  Center  1 

Code  5042 

Crane,  Indiana  47522 
Attn:  Dr.  B.  Douda 

Rohm  and  Haas  Company  1 

723-A  Arcadia  Circle 

Hunsville,  Alabama  35801 
Attn:  Dr.  H.  Shuey 

Strategic  Systems  Project  Office  1 
Dept,  of  the  Navy 
Room  901 

Washington,  DC  20376 
Attn:  Dr.  J.  F.  Kincaid 

Strategic  Systems  Project  Office  2 
Oept.  of  the  Navy 
Room  1048 

Washington,  DC  20376 
Attn:  Mr.  E.  L.  Throckmorton 
Mr.  R.  Kinert 

Thiokol  Chemical  Corp.  Brigham  1 
City 

Wasatch  Division 
Brigham  City,  UT  84302 
Attn:  Dr.  G.  Thompson 
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USA  ARRADCOM  1 

ORDAR-LCE 

Dover,  NJ  07801 

Attn:  Dr.  R.  F.  Walker 

USA  ARRADCOM  1 

DRDAR-LCE 

Dover,  NJ  07801 

Attn:  Dr.  N.  Slagg 

U.S.  Army  Research  Office  1 

Chemistry  Division 

P.O.  Box  12211 

Research  Triangle  Park,  NC 

27709 


Univ.  of  Massachusetts 
Department  of  Chemistry 
Amherst,  MA  03003 
Attn:  Professor  P.  Li  Ilya 

Office  of  Naval  Research 
800  N.-  Quincy  St. 

Arlington,  VA  22217 
Attn:  Dr.  G.  Neece  Code  472 

SRI  International 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 
Attn:  Dr.  D.L.  Ross 


Fluorochem,  Inc.  1 

680  South  Ayon  Ave. 

Azusa,  CA  91702 
Attn:  Dr.  Kurt  Baum 

Univ.  of  Illinois  at  Chicago  Circle  1 
Dept,  of  Chemistry 
Box  4348 

Chicago,  Illinois  60680 
Attn:  Professor  J.H.  Boyer 

University  of  Massachusetts  1 

Department  of  Chemistry 
Amherst,  MA  03003 
Attn:  Professor  J.C.  Chien 

Rockwell  International  1 

Rocketdyne  Division 

6633  Canoga  Avenue 

Canoga  Park,  CA  91304 

Attn:  Dr.  Karl  0.  Christe 

Rockwell  International  1 

Rocke.dyne  Division 
663?  Canoga  Ave. 

Canoga  Park,  CA  91304 
Attn:  Dr.  M.B.  Frankel 


f'olvsciences,  Inc.  1 

Valiev  Industrial  Park 
Warrington,  PA  18976 
Atrn:  Dr.  B.  David  Halpern 

The  Johns  Hopkins  University  1 

Department  of  Chemistry 
Baltimore,  MD  21218 
Attn:  Dr.  Joyce  J.  Kaufman 
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